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Introduction How can the spin 1/2 object known as the 
proton have a non-spherical shape? Why would a physi- 
cist even think of such a concept? Can a non-sphericity 
(or pretzelocity) be measured or computed? This note is 
concerned with such questions. 

The notion that the proton might not have a spherical 
shape has its impetus in the discovery that the spins of 
quarks and anti-quarks account for only about 30% of the 
total angular momentum Many experiments have 
sought the origins of the remainder, expected to arise 
from quark and gluon orbital angular momentum or from 
pairs of strange quarks. 

ij This article is concerned with the relation between the 
quark orbital angular momentum and the non-spherical 
shape of the proton. A number of concerns arise imme- 
diately. As a particle of spin 1/2, the proton can have no 
quadrupole moment, according to the Wigner-Eckart the- 
orem. In elastic electron-proton scattering experiments, 
the effects of relativity cause the initial and final wave 
functions to differ because their momenta differ. For 
example, one thinks of a particle in relativistic motion 
having a pancake shape because of the effects of Lorentz 
contraction. Such an effect is not a manifestation of the 
intrinsic proton shape. 

The presence of significant orbital angular momen- 
tum can only lead to a non-spherical shape if such can 
be defined by an appropriate operator. We used the pro- 
ton model of Ref. to show 0] that the rest-frame 
ground-state matrix elements of spin-dependent density 
operators reveal a host of non-spherical shapes. 

Experimental genesis The electromagnetic current 
matrix element can be written in terms of the Dirac 
F\{Q 2 ) and Pauli F 2 (Q 2 ) form factors, Q 2 is the nega- 
tive of the square of the space-like four-momentum trans- 
fer. These form factors are probability amplitudes that 
the proton can absorb a squared four momentum trans- 
fer Q 2 and still remain a proton. Two exist because 
the rapidly moving quarks within the proton carry both 
charge and magnetization densities. For Q 2 = the form 
factors F\ and kF 2 are the charge and the anomalous 
magnetic moment n in units e and e/2Mj^, and the mag- 
netic moment // = ! + «. The Sachs form factors are 



G E 



kF 2 ,G m = F 1 + kF 2 



In the non-relativistic quark model, Ge and Gm are 
Fourier transforms of the ground state matrix elements 
of the quark charge {J2i=i 3 e i$( r ~ r i) and magnetiza- 
tion J2i=i,3 2m~^( r ~~ r i) density operators. Thus, non- 
relativistically, one expects that G E (Q 2 )/ ^G M (Q 2 ) = 1- 
Interestingly, the opposite highly relativistic limit, in 
which dimensional counting applies, predicted B (us- 
ing the notion of helicity conservation in the inter- 
actions between photons and massless fermions) that 
limg2_ >00 QF2/F1 = m q /Q 2 , where m q is the small mass 
of a down or up quark. This is equivalent to the non- 
relativistic expectation. Thus theoretical expectations 
(and early data) were that the ratio of the Sachs form fac- 
tors would be constant. These expectations were dramat- 
ically thrown aside with the discovery that Ge/Gm falls 
rapidly with increasing values of Q 2 and that QF 2 /F\ 
is approximately constant 0], Q]. See FigJTj which also 
displays the results of our 1995 theory (3,0]. 

An explanation 0] of how the model of 0] describes 
the data showed that the constant ratio QF 2 / F\ emerges 
from the model's relativistic aspects. For the proton 
wave function, only the component in which the first two 
quarks have a vanishing total angular momentum enters 
in computing the electromagnetic form factors. Then the 
angular momentum of the proton S is governed by that 
of the third quark. The relevant Dirac spinor is: 



u(K,S) 



1 



E(K) + m q 



(E(K) + m g )\S) 
a-K\S) 



(1) 



with E(K) = (K 2 + to 2 .) 1 / 2 . The magnetic quantum 
number of the proton is denoted by S, and the lower 
component contains a term a ■ K that allows the quark 
to have a spin opposite to that of the proton's total an- 
gular momentum. The vector K reveals the presence of 
the quark orbital angular momentum: the struck quark 
may carry a spin that is opposite to that of the proton. 
Consequently nucleon helicity 0] is not conserved 10,11]. 



Spin-dependent density operators We interpret or- 
bital angular momentum in terms of the shapes of the 
proton by these are exhibited through the rest-frame 



FIG. 1: The ratio QF 2 /F 1 . The curves from the 1995 theory c 
0| for the ratio are labeled by the value of a model paramete: 
The data are from Q and [f| . Figure reprinted with permissio 
from Q ■ Copyright 2002 by the American Physical Society. 



ground-state matrix elements of spin-dependent densit 
operators 0]. The usual quantum mechanical densit 
operator is p(r) = J2i$( r ~ r i)> where r$ is the positio 
operator of the i'th particle; but for particles of spin 1/2 
one can measure the combined probability that particle 
is at a given position r and has a spin in an arbitrary, 
fixed direction specified by a unit vector n. The resulting 
spin-dependent density SDD operator is 



p(r,n) =5^<5(r- r i)~( l + a i • n )- 



(2) 



To understand the connection between the spin- 
dependent density and orbital angular momentum, con- 
sider a first example of a single charged particle mov- 
ing in a fixed, rotationally-invariant potential in an en- 
ergy eigenstate |*i,i,i/2,s) of quantum numbers: I = 
l,j = 1/2, polarized in the direction s and radial wave 
function R(r p ). The wave function can be written as 
(r p |*l,l,l/2,s) = R(r p )cr • r p \s). The ordinary density. 
p{r) = (*i,i,i/ 2 , s |5(r - r p ) 1*1,1,1/2,5) = R 2 {r), a spher- 
ically symmetric result because the effects of the Pauli 
spin operator square to unity. But the matrix element of 
the SDD is more interesting: 



p(r,n) 



R 2 {r) 



(s\(T ■ f (1 + a ■ n)<r • f |s). 



(3) 



The magnetic quantum defines an axis, s and the direc- 
tion of vectors can be represented in terms of this axis: 
s ■ f = cos#. Suppose n is either parallel or anti-parallel 
to the direction of the proton angular momentum vector 
s. Then p(r, n = s) = R 2 (r) cos 2 6, p(r,n = — s) = 
R 2 (r) sin 2 6, and the non-spherical shape is exhibited. 
The average of these two cases is a spherical shape. 





FIG. 2: (Color online) Shapes of the proton. S is in the 
vertical direction. Left column quark spin parallel to nucleon 
spin. Right column : quark spin anti-parallel to nucleon spin. 
The value of K increases from to 1 to 4 GeV/c. Figure 
reprinted with permission from d. Copyright 2003 by the 
American Physical Society. 



Another useful example is that of the Dirac four- 
component spinor electron wave function of the hydrogen 
atom ground state, with relative size of the lower com- 
ponent governed by the fine structure constant, a. The 
expectation value of the spin-dependent density opera- 
tor, computed using Dirac matrices, with a = 7°7 • n, 
is p(r, h = s) oc [l + a 2 /4cos 2 9] ~ 1 + 10~ 5 cos 2 with 
n = s. For n = — s, p(r, n = — s) = a 2 sin 2 6/4. Relativ- 
ity as manifest by lower components of the Dirac wave 
function causes the hydrogen atom to be slightly, but 
definitely, non-spherical! 

The notion of the SDD can be extended. In condensed 
matter applications [l3| neutrons interact with atomic 
electrons, and only the (electronic) spin-dependent term 
of Eq. ([2|) is used. For quark systems, the densities could 
be weighted by the charge or flavor of the quarks, or use 
other operators. In particular, we use 121 ] 



/W,n 



n 



(4) 



which is more experimentally accessible. 

Now turn to the proton. Its wave function is specified 
in momentum space, so we define [3] a charge-weighted 



3 



SDD operator (the probability that a quark has a mo- 
mentum K and spin direction n): 



p(K,n) 



d 3 £ 



^(£)-(7°+7-n 7 5>K0), (5) 



where Q/e is the quark charge operator in units of the 
proton charge. The quark field operators are evaluated 
at equal time, £° = 0. For the case of /5r the term 7 • n 
is replaced by 7°7 • n. For a spin-polarized (in the S 
direction) proton at rest |^s); the matrix elements of 
p, /5r (SDDs) are given by 

p(K, n, S) = A + Bn • S + C(n • K S • K - -n • S) 

3 



p R (K,n,S) = A R + B R n-S 
+C R (n-kS-K-in-S), 



(6) 



where A,B,--- are scalar coefficients. These forms rep- 
resent the most general rest frame shape of the proton, 
if parity and rotational invariance are upheld 0] • 

We display the shapes of p(K, n, S) 0] in Fig. 2 for 
the cases of quark spin parallel and anti-parallel to the 
polarization direction of the proton S. The shape for a 
given value of K is determined by the ratio of the upper 
to lower components of the quark Dirac spinor Eq. ([1]). 
The relatively large value of the ratio implies considerable 
non-sphericity and a sharp contrast between the proton 
and hydrogen atom. As the value of K increases from to 
4 GeV/c the shape varies from that of a sphere to that of 
a peanut, if n || S. The torus or bagel shape is obtained if 
— n || S. Taking nlS leads to some very unusual shapes 
shown in Fig. 3. Using the given model [5(, one may also 
obtain in coordinate space SDDs. Possible shapes include 
a pretzel form [jj. 

Any wave function yielding a non-zero value of the 
coefficient C(K 2 ) or Cr(K 2 ) represents a system of a 



non-spherical shape. If the relativistic constituent quark 
model of 0] is used, the extra 7 changes the sign of the 
lower component of the wave function, causing Cr = —C. 
Thus either Cr or C can be used to infer information 
about the possible shapes of the nucleoli. Measuring ei- 
ther would require controlling the three different vectors 
n, S and K. 

A specific aspect of p(K, n) is easily related to com- 
pleted experiments because / d?Kp(K., n) is a local op- 
erator. Its matrix element is a linear combination of the 
charge, integrals of spin-dependent structure functions 
Aq (quark contribution to the proton total angular mo- 
mentum), and gA that can be determined from previous 
measurements. We find 



d 3 K{N\p(K, n = S, S) - p(K, n = -S, S)\N) 



0.68, 



(7) 



in which Aq = 0.3 [H and gA = 1.26. The model we use 
gives 0.74 for the above quantity, indicating that shapes 
discussed here may not be unrealistic. 

Measuring the Non-Spherical Shape of the Nu- 
cleoli Can non-spherical shapes be measured? While 
measurements of the matrix elements of the non- 
relativistic spin-density operator [2] reveal highly non- 
spherical densities, finding the non-spherical nature of 
the proton has remained a challenge. Here we explain 
how matrix elements of the spin-dependent density may 
be measured using their close connection with transverse 
momentum dependent parton densities. 

The densities of Eq. © require that the system be 
probed with identical initial and final states. But this 
condition also enters in measurements of both ordinary 
and transverse-momentum-dependent TMD parton dis- 
tributions. The latter lj] are: 



d 2 £,T jK-i 



2 (2vr) J 



(P,S\W)TC$,bn-)1>{JZ)\P,S) 



(8) 



r 



where the specific path n~ is that of Appendix B of [lj]. operator T can be any Dirac operator, e.g. T = icr 
The functions $' r l depend on the fractional momentum \/27°7 l 7 5 , related to hi T , which causes the non-spherical 
x = K + /P + , Ky and on the hadron momentum P. The nature of /)r. 



FIG. 3: Shapes of the proton with n • s = 0. Left column, 
n points (out of page), central: n points sideways, right n is 
out of the page at a 45° angel. The momentum K increases 
from 1 to 4 GeV/c. Figure reprinted with permission from 
Copyright 2003 by the American Physical Society. 



It is therefore tempting to try to associate an SDD 
such as that of Eq. ([5|) with TMDs, but one difference 
is essential. Parton density operators Eq. (JSj) depend on 
quark-field operators defined at a fixed light cone time 
^ + = £ 3 + £° = while our SDD is an equal-time, £° = 0, 
correlation function. However, a relation between the two 
sets of operators is obtained [L2J by integrating the TMD 
over all values of x setting £~ to zero, and integrating 
Eq. (0) over all values of K z so that £ 3 = 0. After integra- 
tion, £ — for both functions. The density operators, 
derived from those of Eq. §6§ are denoted by adding a T 
to the subscript. Thus /3Rr(KT,n) = /f^, dK z p^(K, n). 
It is therefore a matter of algebra to show that 

p RT (K T ,n T , S T )/M = h{K T ) + hx(K T )n ■ S T 
(n T • K T S T • K T - ±Kfn • S T ] 



+ 



M 2 



-hi T {K 2 T ), 



(9) 



in the rest frame, where a tilde is placed over each TMD 
parton distribution to denote an x-integrated function. 
Finding that non-zero value of h^ T 7^ would demon- 
strate that the proton is not spherical. 

The term h^ T causes distinctive experimental signa- 
tures in semi-inclusive leptoproduction hadron produc- 
tion experiments 15, Q3]. If the target is polarized in 



a direction St transverse to the lepton scattering plane, 
the cross section acquires a term proportional to cos (3<p l h ) 
where (j) l h is the angle between the hadron production 
plane (defined by the momenta of the incoming virtual 
photon and the outgoing hadron) and the lepton scatter- 
ing plane. A similar effect occurs in electroweak semi- 



inclusive deep inelastic leptoproduction 17J. In each of 
these cases, the momentum of the virtual photon and its 
vector nature provide the analogue of the vector n needed 



FIG. 4: Transverse shapes of the proton; 

V^/Sri^Kt, n)//i(if|i). The horizontal axis is the the 
direction of St and n = St, 4> n = it. The shapes vary from 
circular to highly deformed as Kt is increased from to 2.0 
GeV in steps of 0.25 GeV. Figure reprinted with permission 
from Copyright 2007 by the American Physical Society. 



to define the spin-dependent density. The hadronic trans- 
verse momentum provides the third, Kt- Another possi- 
bility occurs in the Drell-Yan reaction pp(]) — ► UX, using 
one transversely polarized proton 

The shapes inherent in Eq. Q are illustrated using the 
spectator model of 19)]. Here <f> is the angle between Kt 
and St and (fr n is the angle between n and St • The trans- 
verse shapes of the proton (assuming a struck u quark) 
are shown in Fig. 4, taking <j) n = ir. This emphasizes 
the non-spherical nature because the first two terms of 
Eq. ([9]) tend to cancel. The shapes of Eq. (|9|) can be 
thought of as projections of the shapes displayed in pre- 
vious figures. 



The model 19[] indicates that the functions f\ , h\ and 
hj- T have very similar x dependence, so that measure- 
ments at values of x for which these functions peak should 
be sufficient to construct the required integrals over x. 

The non-spherical nature of the nucleon shape is de- 
termined by the non- vanishing of the TMD hj~ T . It is 
very exciting that experiments planned at Jefferson Lab- 
oratory aim to specifically measure h^ T 20] and therefore 
determine whether or not the proton is round. 

Connection with lattice QCD The non-spherical 
shape of the nucleon can be established in lattice QCD 
by computing the lattice version of the angular integral 
of the matrix element: 

F r (r) = [ dfy 20 (r)(^s|VKr)(7 + ^ • n 75 )V'(0)|^ s ) 



where T = 1 or 7 and the link operator is not displayed. 



5 



A non-zero value of F(r) for any value of r would imme- 
diately tell us that the proton does not have a spherical 
shape. Matrix elements of ^(^^(O) have been evalu- 
ated for the case when the separation is one or two links. 
Thus the relevant information is available. Preliminary 
results for -Fr( r ) exist only for separations of one-link, 
and current statistics are not high [21[. Another possi- 
bility, closely related to finding h^ T , would be to take 
the spatial component of r to be perpendicular to s and 
integrate over the transverse directions. I hope that the 
lattice QCD community will find it of sufficient interest 
to warrant the effort of a detailed, high-statistics calcu- 
lation. 

Summary The nature of the proton wave function can 
be elucidated by studying the matrix elements of a gener- 
alized density operator. Spin-dependent quark densities 
SDD are defined as matrix elements of density operators 
in proton states of definite spin-polarization, and shown 
to have an infinite variety of non-spherical shapes. For 



high momentum quarks with spin parallel to that of the 
proton, the shape resembles that of a peanut, but for 
quarks with anti-parallel spin the shape is that of a bagel. 
The matrix elements of the SDDs are closely related to 
specific transverse momentum dependent TMD parton 
distributions accessible in the angular dependence of the 
semi-inclusive processes ep — > enX and the Drell-Yan re- 
action pp — > UX. New measurements or analyses would 
allow the direct exhibition of the non-spherical nature 
of the proton. The TMDs can be computed using lattice 
QCD so that the non-spherical shapes could be measured 
experimentally and computed using fundamental theory. 
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